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We study magnetically ordered phases and their phase boundaries in the J\ — J2 — J3 Heisenberg 
models on the honeycomb lattice using series expansions around Neel and different colinear and 
non-colinear magnetic states. An Ising anisotropy (A = J±/J z 7^ 1) is introduced and ground state 
energy and magnetization order parameter are calculated as a power seies expansion in A. Series 
extrapolation methods are used to study properties for the Heisenberg model (A = 1). We find that 
at large J3 (> 0.6) there is a first order transition between Neel and columnar states, in agreement 
with the classical answer. For Jz — 0, we find that the Neel phase extends beyond the region 
of classical stability. We also find that spiral phases are stabilized over large parameter regions, 
although their spiral angles can be substantially renormalized with respect to the classical values. 
Our study also shows a magnetically disordered region at intermedaite J2/J1 and J3/J1 values. 

PACS numbers: 74.70.-b,75.10. Jm,75.40.Gb,75.30.Ds 



INTRODUCTION 



The search for quantum spin-liquid phases in realistic 
spin and electronic models and real materials remains a 
very active area of research. In recent years a number 
of experimental materials have been synthesized, whose 
behavior has been very promising from the point of view 
of discovering such a phase. These include many kagome 
and triangular lattice based frustrated magnets. [l|-|4| In 
several cases, the frustration parameter, [5[ defined as the 
ratio of the Curie- Weiss temperature to any magnetic 
ordering temperature exceeds 1000. Most of these exper- 
imental systems show gapless spin excitations, though 
the role of impurities and weak anisotropics in creating 
gapless spin excitations has not been ruled out. 

From a theoretical point of view, frustration is es- 
sential for obtaining a spin-liquid phase in dimensions 
greater than one. Unfrustrated two-dimensional spin 
models show robust Neel order. However, both spin- 
frustration arising from competing exchange interactions 
and frustration arising from the itineracy of the electrons 
can lead to a spin- liquid phase. In recent years several 
realistic models have emerged as candidates for quan- 
tum spin-liquid ground states. One is the Heisenberg 
model on a kagome lattice, where extensive Density Ma- 
trix Renormalization Group work by Yan et alQ presents 
strong evidence for a spin-liquid phase with a not-too- 
small gap to singlet and triplet excitations. Another 
is the triangular-lattice antifcrromagnet with Heisenberg 
and 4-spin ring exchanges. [2IIH A third, somewhat unex- 
pected case, is the half-filled honeycomb lattice Hubbard 
model, Q which would lead to an unfrustrated Heisenberg 
model at large U. It shows a spin-liquid phase at inter- 
mediate values of U, sandwiched between a Neel phase 
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FIG. 1: Exchange constants of the honeycomb-lattice Heisen- 
berg model 



at large U and a semi-metal phase at small t/.[l(J This 
spin-liquid phase was also found to be gapped with no 
signs of any broken symmetry. 

The question of whether this latter phase can be re- 
alized in a frustrated spin model on the honeycomb lat- 
tice, arising from higher order charge exchange processes 
has attracted considerable interest. [11| Several recent pa- 
pers have investigated this question using exact diagonal- 
ization, Schwinger-Boson and spin-wave theory, coupled 
cluster methods, and variational wavefunctions. They 
have come to conflicting conclusions about the phase di- 
agram and the possible existence of a spin-liquid phase in 
this model. [12h18| Experimental realizations of spin- half 
honeycomb lattice Heisenberg model have been discussed 
recently by Tsirlin et al. 19] 

Here we study the phase diagram of frustrated Heisen- 
berg model on the honeycomb lattice using series 
expansions. [13, 0] We will confine our study to all ex- 
changes being positive (Ji, J2, J3 > 0) (See Fig. 1). The 
nearest-neighbor honeycomb lattice Heisenberg model, 
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FIG. 2: Classical phase diagram of the honeycomb lattice 
Heisenberg model with J\ = 1. 



dimer state is called a staggered dimer state. However, 
we found that the staggered dimer phase has rather high 
energy and thus is not stabilized for J 3 = 0. While our 
study can not rule out such a phase for J 3 = 0.3 and 
intermediate J2, the energies again suggest that the stag- 
gered dimer state has too high an energy. This should be 
contrasted with studies of the square-lattice J± — J2 mod- 
els, where energies of dimer-phases match smoothly with 
the magnetically ordered phases at intermediate J2 / J\ 
values. |30j It is possible that some other dimer or plaque- 
tte phase is stable here.[18[ Investigation of larger unit 
cell Dimer/Plaquette phases is left for future work. 



MODELS AND SERIES EXPANSIONS 



with antiferromagnetic exchange J\ (which we set to 
unity), is well known to have a Neel ordered ground 
state. [2UH23l |25| . |26 | Adding a frustrating second neigh- 
bor antiferromagnetic interaction J2 and a third neigh- 
bor antiferromagnetic interaction J3 leads to a complex 
phase diagram even at a classical level 2J, |26| that in- 
cludes Neel, columnar and different spiral phases (See 
Fig. 2). We have carried out series expansions around 
Neel, columnar and spiral phases. For large J3 (> 0.6) 
we find a direct first order transition between Neel and 
columnar phases as a function of J2. Thus, the most in- 
teresting part of the phase diagram is for J3 < 0.6. We 
have studied this phase diagram along constant J3 lines 
as a function of J2 as well as along the contour J3 = J2 . 
We find that at small J3 the Neel phase widens due to 
quantum fluctuations, in agreement with recent exact di- 
agonalization study. [18| We also find that spiral phases 
are stabilized for intermediate and large J 2 . The spi- 
ral angles are strongly renormalized with respect to their 
classical values. However, they do not lock into commen- 
surate colinear structures as suggested in recent exact 
diagonalization study. [3] Rather, our studies favor an- 
gles closer to 90 degrees, which implies spins pointing at 
right angles. 

Our study also shows that a part of the phase diagram 
at small and intermediate J3 and J2 has no magnetic long 
range order. The convergence of our series expansion is 
best for colinear phases, when they have clearly non-zero 
order parameters. We are unable to accurately obtain 
the region of stability of magnetically disordered phases, 
but it is clear from the general analysis that such a phase 
must exist. This is also clearly seen by going along J 2 — 
J3, where such a phase exists near the highly frustrated 
point J 2 — J3 = 0.5. 13] It is possible that a strip of this 
phase extends all the way down to J 3 = 0. 

In order to gain some insight into the nature of this 
phase, we have also carried out series expansions around 
a dimer state, which was considered as a possible ground 
state in Ref. [l2|. Following Moessner et al,[29j such a 



We study the honeycomb-lattice Heisenberg model 
with Hamiltonian 

H = Ji & ■ Sj + ,h J2 Si ■ Sk + J 3 J2^- Si- (1) 
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The exchanges within an elementary hexagon are shown 
in Fig. 1. Without loss of generality, we take J\ = 1. The 
classical phase diagram for J 2 > 0, J3 > is sketched in 
Fig. 2. The spin structure of the classical phases are 
shown in Fig. 3. It consists of the Neel phase, a colum- 
nar phases, which is also colinear, and two planar spiral 
phases. In the spiral-1 phase the pitch vector Q is perpen- 
dicular to one of the bond directions, and is characterized 
by a single angle 9. Spins twist away from the Neel state 
by this angle as one goes from one row to the next. In 
the spiral-II phase the pitch vector is parallel to one of 
the bond directions. It is characterized by two angles, 
9 and 4>. Within a unit cell, the spins deviate from the 
Neel state by an angle 4>. In addition, there is a twist 
by an angle theta as one goes from one unit cell to next. 
The spiral II phase is stabilized by negative J3 but be- 
comes degenerate with the spiral I phase for J3 = and 
1/6 < J 2 < 1/2. Also, at J 2 = 1/2, the spiral-II phase 
locks into another colinear phase with (j) = tt, 9 = n. 

To carry out an Ising expansion around a colinear state 
(Neel or columnar), we introduce an Ising anisotropy in 
all the exchange interactions, 



Si ■ Sj — S*Sj 



\{S*S* + SfS]). 



(2) 



We calculate series expansions in the variable A for the 
ground state energy, per spin, eo and the sublattice mag- 
netization m. 

The reader unfamiliar with series expansion methods 
can find technical details in refs. 25,26 and we shall not 
elaborate here. However it is worth mentioning that it 
is convenient, for both colinear and non-colinear states, 
to rotate the spin axes on different sites to give an un- 
perturbed ground state in which all spins are aligned. 
This procedure has been described in ref. [32j. This trans- 
formation has the effect of complicating the form of the 
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FIG. 3: Classical phases of the honeycomb-lattice Heisenberg 
model with positive Ji, J2 and J3. The phases are (a) Neel, 
(b) columnar, (c) spiral-I and (d) spiral-II. 





FIG. 5: Ground state energy estimates of the Neel phase 
(J2 < 0.5) and the columnar phase (J2 > 0.5) are shown 
as a function of J2 for different values of J3. For J3 = 0.6 
and J3 = 1.0 linear fits to the data points are also shown 
to illustrate the first order transition point between the two 
phases. 



FIG. 4: Staggered dimer phase on the honeycomb lattice 



Hamiltonian but avoids the need to identify sublattice 
types in the cluster data. 

To derive series to a given order n requires enumera- 
tion of clusters with up to n sites. This number grows 
rapidly with n, and this is the major factor limiting the 
length of series. For example, for the 8th order series for 
the spiral I phase for the full J1-J2-J3 model a total of 
1,083,315 distinct clusters occur. For the J1-J2 model 
an additional term can be obtained, requiring a total of 
1,172,204 clusters. 

The spiral II phase is only considered for J3 = (it is 
not stabilized classically for positive J3). It has 7 bond- 
types and took 108,453 distinct clusters to calculate the 
series to 8th order. For both the spiral phases, we allow 
angles ( such as 8, (f>) different from classical values and 
take the angles that minimize the ground state energy. 

In addition, we have also developed series expansion for 
the ground state energy of the dimer configuration shown 
in Fig. 4. The center of the dimers form a triangular- 
lattice with three different types of inter-dimer interac- 
tions with Ji, J r 2 and J3 all non-zero. To carry out the 
calculation to 8-th order required 33826 linked clusters. 

The series data is too lengthy to reproduce here, but 
can be provided on request. 



SERIES ANALYSIS AND RESULTS 

We will restrict attention to a few contours in the J2 — 
J3 plane. We will consider J 3 = 0, J3 = 0.3, J3 = 0.6 
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FIG. 6: Magnetization for Neel phase (J2 < 0.5) and colum- 
nar phase ( J2 > 0.5) as a function of J2 for different values of 
J3 . The arrows indicate the critical J2 values for the classical 
model, where the Neel order parameter vanishes, for J3 = 0.0 
and J3 = 0.3. 



and J3 = 1 , and J3 = J2 ■ For large J 3 Neel and columnar 
states are considered, whereas for small J3 Neel and spiral 
phases are considered. 

The series are analyzed by Pade as well as d-log Pade 
approximants. These approximants are constructed in 
the variable A as well as in the variable 8 = 2A — A 2 . The 
latter (3^] allows one to eliminate a square-root singular- 
ity due to the gapless spin- waves at A = 1. 

In Fig. 5, we show the ground state energy of the 
Neel and columnar phases for different values of J3 . The 
corresponding magnetizations are shown in Fig. 6. For 
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FIG. 7: Order parameter for the Neel and columnar phases 
along J3 — J2. 
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FIG. 8: Ground state energy for Neel, staggered dimer, spiral- 
I and spiral-II phases for J3 = 0. 



Phase 


J2, J3 


del, 4>cl 


9 r , 4>r 


Spiral-I 


1.00, 0.0 


104.5 


90.0 


Spiral-I 


0.70, 0.0 


98.2 


90.0 


Spiral-I 


0.50, 0.0 


90.0 


85.0 


Spiral-I 


0.70, 0.3 


104.5 


95.0 


Spiral-I 


0.60, 0.3 


99.6 


80.0 


Spiral-II 


0.40, 0.0 


149,125 


150,95 


Spiral-II 


0.25, 0.0 


104.5,75.5 


120,80 



J =0.3 



TABLE I: Selected examples of classical spiral angles ( 9 c i 
for spiral-I phase and 9 c i and <j> c i for spiral-II phase) and the 
corresponding estimated renormalized spiral angles (8 r and 
4>r) in degrees. 



J 3 = 1.0, there is a clear first order phase transition be- 
tween the Neel and columnar phases around J2 > 0.5. An 
estimate of the transition point can be obtained by the 
linear fits to the energy data (also shown). For J3 = 0.6, 
the transition from the two sides is close to becoming con- 
tinuous. The order parameters for both phases approach 
zero as J2 approaches 0.5. Very near the transition, the 
series analysis is not reliable enough to tell if this is a 
continuous transition, or a first order transition, or there 
is a small intermediate region with no magnetic order, 
although a direct transition between the phases will nec- 
essarily be first order, as seen by the linear fits to the 
ground state energy. 

The sharp downturn in the columnar-state magnetiza- 
tion for larger J2 is a clear signature of an impending 
transition to the spiral I phase, which classically would 
occur at J2 = 0.7. 

In Fig. 7, we show the magnetization for the Neel and 
columnar phases along the contour J2 = J3 in the param- 
eter space. We can see that both order parameters go to 
zero before the highly frustrated point J2 = J3 = 0.5 is 
reached and there is an intermediate phase with no mag- 
netic order. This agrees with previous theoretical study 








FIG. 9: Ground state energy for Neel, staggered dimer and 
spiral-I phases for J3 = 0.3. 



by Cabra et al. [13| 

At small J3 (for example J3 = 0) the region of Neel or- 
der is increased with respect to the classical answer, but 
that is no longer true for J3 = 0.3. For J3 = 0, we esti- 
mate that Neel order extends upto J2 = 0.20 compared 
to the classical value of 1/6. In contrast, for J3 = 0.3, 
the magnetization vanishes near J2 = 0.3 compared to 
the classical value of 11/30 (See Fig. 6). We do not find 
any clear evidence for a magnetically disordered phase at 
J3 = 0. Fig. 8 shows the ground state energy for Neel, 
spiral-I, spiral-II and staggered dimer phases for J3 = 0. 
In the spiral phases, we consider a range of spiral an- 
gles and pick that value which minimizes the energy. In 
general, these energy functions are quite shallow, so that 
the renormalized spiral angles shown in Table-I should be 
considred as approximate. While the convergence is not 
excellent at intermediate J%j J\ values, our analysis sug- 
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gests that one transitions from Neel to spiral-II to spiral-I 
phase. The staggered dimer phase is clearly not stabi- 
lized here in contrast to recent suggestions from bond- 
operator formalism[l2j based studies. In addition, the 
spiral-II phase also does not lock into a colinear configu- 
ration in contrast to the conclusions from the finite-size 
study. [lH If anything, our estimates suggest spiral angles 
close to 90 degrees, which implies some spins pointing at 
right angles with respect to their neighbors. 

In Fig. 9, we show the ground state energy for Neel, 
spiral-I, and staggered dimer configurations along J 3 = 
0.3, as a function of Ji- At intermediate Ji values the 
convergence of both Neel and spiral-I phases becomes 
poor. The energy of the staggered dimer phase appears 
well behaved. In this case, there is a much stronger case 
for an intermediate phase with no magnetic order. How- 
ever, the staggered dimer phase has a relatively high en- 
ergy and it is not clear that it becomes the ground state of 
the system in this region. Since the ground state energy 
must vary continuously, it seems more likely that some 
other phase not considered by us becomes the ground 
state in this region. At large J2, the energy of spiral 
and staggered dimer phases become nearly degenerate. 
Our results are consistent with the suggestion from the 
exact diagonalization study [3] that for some parameter 
ranges, the staggered dimer state may be difficult to dis- 
tinguish from magnetically ordered states. 



This difference may be because periodic boundary con- 
ditions on finite systems may disfavor incommensurate 
phases. Another difference with respect to the Bond Or- 
der Mean-Field Theory study [IH is that we do not find 
evidence for a staggered dimer phase in the magnetically 
disordered region at small J3. 

An important question is whether this magnetically 
disordered phase in the frustrated Heisenberg models is 
related to the one found by Meng et al[§| in the Hubbard 
model at intermediate 17/ 1, and whether it is a true spin- 
liquid. Recent work by Yang and Schmidt [HI found that 
for U/t values where Meng et al obtained a phase tran- 
sition to a spin-liquid phase J2/J1 remains very small 
(of order 0.06), with J3/J1 even smaller. This means 
that just these additional interactions can not drive the 
transition, as our study finds that these parameters are 
well within the Neel phase. It is possible that larger 
ring exchanges play a role in bringing about the transi- 
tion. Nevertheless, the disordered phase in the Heisen- 
berg model may be connected to the one found in the 
Hubbard model. This question as well as the possibility 
of other types of dimer/plaquette phases deserves further 
attention. 

This work is supported in part by NSF grant number 
DMR-1004231. 



DISCUSSIONS AND CONCLUSIONS 

In this paper, we have used series expansion meth- 
ods to study the ground state phase diagram of the 
J\ — J2 — J3 antiferromagnetic Heisenberg models on the 
honeycomb lattice. We have determind the stability of 
Neel, columnar and various spiral phases and calculated 
their properties. In agreement with previous studies, we 
find that the region of stability of Neel phase increases 
with quantum fluctuations for J3 = 0. However, that is 
not the case for larger J3. There clearly exists a param- 
eter region at intermediate J2 and J3 with no magnetic 
order. This region is most clearly seen near the highly 
frustrated point J3 = J2 = 0.5 and possibly forms a strip 
extending all the way near to J3 = 0. We found that 
the staggered dimer order is not favored for the J\ — J2 
model. 

The phase diagram has some resemblence to the 
square-lattice Heisenberg model with frustrated antifer- 
romagnetic exchange constants. In all these systems, it 
remains difficult to determine the nature of the magneti- 
cally disordered phases. Our results are generally in good 
agreement with recent exact diagonalization study of Al- 
buquerque et al.[3] One point of difference is that we 
do not find that a third colinear phase (that is the spi- 
ral II with <j) — 7r, 9 — 7r) is favored over general spi- 
ral phases over any extended part of the phase diagram. 
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